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Технології високопродуктивних обчислень 
для дослідження та розв'язування задач 
розрахунку міцності конструкцій 


Розглядаються технології високопродуктивних обчислень для дослідження і розв'язування задач 
(статичних лінійних і нелінійних та динамічних лінійних) розрахунку міцності конструкцій. Описано 
паралельні алгоритми формування, дослідження та розв'язування дискретних задач, які одержуються в 
результаті дискретизації методом скінченних елементів вихідних задач, а також реалізація запропонованих 
технологій в програмному комплексі розв'язування задач розрахунку міцності будівельних конструкцій. 


Вступ 


При проектуванні об'єктів різного призначення часто виникає необхідність в 
проведенні розрахунків міцності конструкції таких об'єктів або окремих її елементів. 
В першу чергу задачі розрахунку міцності конструкцій виникають у промисловому 
та цивільному будівництві - при розрахунках споруд в цілому або їх окремих кон- 
структивних елементів, а також в різних галузях машинобудування - суднобуду- 
ванні, авіабудуванні, ракетобудуванні, моторобудуванні тощо. 

Зростаючі вимоги до якості проектних рішень, застосування нових конструк- 
тивних матеріалів викликає необхідність у розв'язуванні якісно нових задач. Зростає 
потреба в нових методах 1 підходах, пов'язаних з побудовою 1 дослідженням корект- 
них комп'ютерних моделей, які адекватно відображають реальну поведінку кон- 
струкцій. Зростають вимоги до достовірності одержуваних комп'ютерних результатів. 
Ці чинники призводять до того, що істотно зростає обсяг оброблюваної інформації. 
Також великий обсяг інформації необхідно обробляти при виконанні розрахунків уні- 
кальних конструкцій. Тому виникає потреба у використанні високопродуктивних об- 
числень, що у свою чергу викликає зростання вимог до обчислювальних ресурсів. 

З іншого боку, на сучасному етапі розвитку обчислювальної техніки відбувається 
перехід від комп'ютерів з одним одноядерним процесором, можливості підвищення 
продуктивності яких вже вичерпано, до комп'ютерів з багатоядерними процесорами, 
яких теж може бути багато. Використання багатоядерних і багатопроцесорних комп'ю- 
терних систем є необхідною передумовою створення технологій високопродуктивних 
обчислень. Але ефективність таких технологій визначається переважно відповідністю 
між програмно-алгоритмічним забезпеченням і технічними засобами багатоядерного і 
багатопроцесорного комп'ютера. Адже найвищої продуктивності можна досягти, ви- 
користовуючи алгоритми і програми з паралельною організацією обчислень. Тому 
виникає потреба як в адаптуванні для багатопроцесорних комп'ютерів раніше створе- 
ного програмного забезпечення комп'ютерів з одним одноядерним процесором, так і в 
створенні якісно нового програмного забезпечення, яке використовує паралельні 
алгоритми дослідження та розв'язування задач. 

Метою цієї роботи є впровадження високопродуктивних обчислень в існуючі 
програмні засоби для комп'ютерів з одним одноядерним процесором. 
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Так, для проведення розрахунків міцності конструкцій створено багато про- 
грамних засобів, як універсального призначення, так і орієнтованих на певну галузь 
(ХАЗТКАМ, МАВС, АМ5У5 тощо). Наприклад, для розрахунку міцності будівель- 
них конструкцій розроблено вітчизняний програмний комплекс (ПК) ЛІРА (|1|. На- 
строївши його інтерфейс для відповідної предметної області, ПК ЛІРА можна вико- 
ристовувати для розрахунків в інших галузях. 

На базі препроцесора 1 постпроцесора ПК ЛІРА та інтелектуальних паралельних 
програм бібліотеки Іпрагіїб (2), |3| розроблено програмний комплекс, який реалізує 
технології високопродуктивних обчислень для дослідження і розв'язування задач роз- 
рахунку міцності будівельних конструкцій на комп'ютерах з паралельною організа- 
цією обчислень, зокрема, на інтелектуальних паралельних комп'ютерах Інпарком |4). 


Математичні постановки задач 


Математично задачі розрахунку міцності конструкцій з використанням прин- 
ципу можливих зміщень можуть бути поставлені у вигляді варіаційних задач: необ- 
хідно знайти вектор-функцію иє (о, яка для будь-якої вектор-функції ує (/о (для 
будь-якого можливого зміщення) задовольняє одній з інтегральних тотожностей 

-- для статичної задачі 


а(и,у) - У); 1) 
- для динамічної задачі 
а(иуу) 1 Б(и"уу) - с(и'уу) - ЦД у), | (о) - З ц'(вугиб»; (2) 
- для задачі на власні коливання 
а(и,у) - 25(и,у), (3) 


де (о - нескінченновимірний функціональний простір можливих переміщень, си- 
метричні білінійні функціонали а(и,у) - Ї єї (») Чи)аЙ , Б(и"м»), с(и?у) пропорційні 
у 
відповідно потенційній, кінетичній енергіям деформації і енергії гальмування, а лі- 
нійний функціонал /(/,у) пропорційний роботі прикладених (зовнішніх) зусиль при 
навантаженні (тут и! З (их, цу, И,) - вектор-функція зміщень, через и "позначено її пер- 
шу похідну за часом, а через и" - другу, є З | Єкк Єуу» Є» Єху» бух Єхх| - ВЕКТОР Дефор- 
мацій, т з Ї ух» Туу» То» Тху» Ту ТухІ - ВеКТОр напруг, Й - об'єм, 5 - поверхня тіла). 
Вектори напруг т і деформацій є зв'язані такою рівністю 


т - Св, (9) 


де С - узагальнена матриця пружності. У випадку нелінійної задачі елементи цієї 
матриці С залежать від деформацій, а отже, і від зміщень, на відміну від випадку 
лінійної задачі, коли ці елементи не залежать від шуканої функції и. 

Теоретичною основою більшості програмних засобів для розрахунку міцності 
конструкцій є метод скінченних елементів (МСЕ), реалізований у формі зміщень. 
Вибір саме цієї форми пояснюється простотою її алгоритмізації і фізичної інтерпре- 
тації, наявністю єдиних методів побудови матриць жорсткості і векторів наванта- 
жень для різних типів скінченних елементів, можливістю врахування довільних гра- 
ничних умов і складної геометрії конструкції, що розраховується. 

Дискретизація варіаційної задачі методом скінченних елементів, який є типовим 
представником проекційних методів, полягає в заміні нескінченновимірного простору 
допустимих функцій (о Його скінченновимірним підпростором ФЛ Вектор-функції з 
підпростору (о можуть бути представлені у вигляді лінійної комбінації базисних 
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вектор-функцій, які задовольняють головним (кінетичним) граничним умовам, 
п 
п 
и"(Х)- Ухеіх) , (5) 
і-«1 


де ф (і - п) - згаданий вище базис Іл звичайно кусково- -поліноміальний. 

б ак і лінійний функціонали з П1)-(3)в скінченновимірному просторі 
ОО у випадку вихідної лінійної задачі можна записати у вигляді відповідно білі- 
нійних і лінійної форм коефіцієнтів х; в представленні (5), які надалі називатимемо 
вузловими нараметрами. 


ауру в у Ах, Бу? зу!Вх, с(и?, у" 2 у'сх, ру! у уть, (6) 
де хіу- вектори вузлових параметрів відповідно функцій и! і у", а елементи 
матриць жорсткості (А), мас а демпфування (С) 1 вектора навантажень (Р) обчис- 
люються за формулами (ії, у - п) 

ду - - а(Ф., Фу), фа бе ун Фу), Су - с(Фь Фу), фіз І; Фі). (7) 
Таким чином, враховуючи в (6) довільність вектора у, отримуємо відповідні 
дискретні задачі |31: 


-- для статичної задачі (1) - систему лінійних алгебраїчних рівнянь (СЛАР) 


Ахаіб; (8) 
-- для динамічної задачі (2) - задачу з початковими умовами 
Вх'"(2) Сх) б Ах() з БО), о х(о)- хо, ху схо (9) 


- для задачі на власні коливання (3) - узагальнену алгебраїчну проблему 
власних значень 
Ах - Х'Вх, (10) 
Але для випадку нелінійної задачі дискретну задачу так отримати не можна. 
Задачу спершу треба лінеаризувати. Як правило, лінеаризація полягає в тому, що 
будується ітераційний процес: вибирається початкове наближення Со для узагаль- 
неної матриці пружності С з (4), тоді (у випадку стаціонарної задачі) отримуємо 
лінійну задачу а ибу") ЩІ; у?), розв'язок якої и! використовується для обчислення 
наступного наближення С! узагальненої матриці пружності. Це наближення вико- 
ристовується для побудови наступної лінійної задачі ази, У бе зу ЗЕ і т.д. доти, доки 
не буде виконана умова збіжності, наприклад, така 


п Й 
і, У | 


«Є. 
Ь Є 
и, 


Таким чином, Б нн процес полягає в знаходженні послідовності функцій 
є (1, 1-1,2,..., кожна з яких при довільній вектор- -функції у"є Ц/ задовольняє 
інтегральній тотожності аа", у «Цру") з білінійним і лінійним функціоналами. 
Кожна з цих лінійних варіаційних задач легко зводиться до СЛАР виду (8) Ах- Б, 
:- 1,2,.... В деяких випадках нелінійність моделюється змінними навантаженнями, 
і тоді ітераційний процес будується так, що від ітерації до ітерації змінюється тільки 
права частина системи (8), а матриця залишається постійною, тобто послідовно 
розв'язуються задачі Ах - Б, І - 1,2, 


Формування і розв'язування дискретних задач 


Враховуючи, що базисні функції підпростору МСЕ ОП! можна вибрати так, щоб 
вони були відмінними від нуля лише на декількох скінченних елементах, структура 
матриць з (6) є в загальному випадку розрідженою і визначається нумерацією вузло- 
вих параметрів, яка у свою чергу залежить від нумерації вузлів скінченноелементної 
сітки. Тобто матриці задач (8) - (10) є стрічковими, профільними, блочно-діагональ- 


«Штучний інтелект» 3'2010 143 


Хіміч О.М., Попов О.В., Полянко В.В. 


ними з обрамленням тощо. Крім того, з точки зору комп'ютерного розв'язування цих 
задач істотне значення має те, що матриці дискретних задач є симетричними, додатно 
означеними або додатно напівозначеними, а також, що порядок цих матриць складає 
0107) - 0(107). 

Формування дискретних задач. Для побудови простору МСЕ (о! вихідна об- 
ласть розбивається на підобласті-елементи. На кожному елементі функції з простору 
МСЕ є поліномами заданого степеня, коефіцієнти яких визначаються через вузлові 
параметри - значення функцій (а у ряді випадків 1 їх похідних) у вузлах елемента 
(точки у вершинах, на сторонах, гранях або усередині елемента). Для вузлових пара- 
метрів і базисних функцій справедливі такі співвідношення 


ІКфі) 7 бу, (11) 
де ді - символ Кронекера, а результатом операції 1/(Фі) є значення компоненти 
вектор-функції Ф; для вузлового параметра /,. 

Важливою перевагою МСЕ, враховуючи (5) 1 (11), є те, що елементи матриць і 
векторів правих частин задач (8) - (10) отримують підсумовуванням відповідних еле- 
ментів відповідних матриць 1 векторів навантажень, побудованих для окремих скінчен- 
них елементів. Така властивість матриць 1 векторів дискретних задач МСЕ дозволяє 
ефективно розпаралелювати процес формування цих матриць і векторів |6|. При цьому 
можливо розпаралелювати тільки обчислення та розподіл між паралельними проце- 
сами відповідно до вимог методу, який використовуватиметься для розв'язування дис- 
кретної задачі, елементів глобальних матриць 1 векторів дискретних задач, використовую- 
чи обчислені раніше матриці і вектори окремих скінченних елементів. Або розпаралелю- 
вати також обчислення цих матриць і векторів окремих скінченних елементів. 

З метою мінімізації кількості арифметичних операцій при розв 'язуванні дис- 
кретних задач, а також необхідної оперативної і зовнішньої пам'яті, враховуючи 
розріджену структуру матриць задач (8)- (10), в більшості випадків проводиться 
переупорядковування невідомих задач. Залежно від критерію мінімізації - ширини 
стрічки матриці, її профілю, загальної кількості арифметичних операцій при трикут- 
ному розвиненні матриці і т.д. - використовуються різні методи переупорядко- 
вування, наприклад, зворотний алгоритм Катхілла-Маккі, фактор дерев, мінімального 
степеня |7|. Конкретних рекомендацій для вибору методу впорядкування дати не 
можна, тому що ефективність того або іншого алгоритму суттєво залежить від пер- 
винної структури конкретної матриці. 

Розв'язування систем лінійних алгебраїчних рівнянь. Для розв'язування 
СЛАР (8) - вихідної або переупорядкованої - на МІМР- комп'ютері можуть бути 
застосовані паралельні алгоритми для дослідження і розв 'язування СЛАР з симетрич- 
ними матрицями, які використовують ІБ грозвинення (або С "рИ-розвинення) ме- 
тоду Холецького |3|, |З|. Залежно від ширини і заповнення стрічки матриці викори- 
стовуються: 

- блочні паралельні алгоритми для вузьких стрічкових матриць і блочно-діаго- 
нальних матриць з обрамленням; 

- одновимірні блочні циклічні паралельні алгоритми для дослідження і роз- 
в'язування СЛАР із стрічковими або профільними матрицями; 

те одновимірний блочний циклічний паралельний алгоритм (ОТ 00-розвинення 
у формі зовнішніх добутків) для дослідження і розв'язування СЛАР з матрицями так 
званої хмарочосної структури (у яких не беруть участь в обчисленнях внутрішні 
елементи профілю матриці, які залишаються нульовими в процесі розвинення матриці). 

Розв'язування задачі з початковими умовами. В напівдискретному МСЕ набли- 
жений розв'язок шукається у вигляді (5), де коефіцієнти х; є функціями часу /. 
В результаті отримуємо систему (9) звичайних диференціальних рівнянь другого по- 
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рядку з початковими умовами, де: х(І), х(, х 


(0). (0) (0) - ( 
х'єб(и"?), хо'еФ(и"). 
В переважній більшості програмних засобів розрахунку міцності конструкцій 
система (9) розв'язується методом розвинення за формами власних коливань. Якщо 
Ж 


2! (2й'В2! 21, К -1,2,...,п) - розв'язок алгебраїчної проблеми власних значень 


1 
Ф вектори з елементами ха), 


п 
(10), то, вважаючи в (9) х(г) - У уц (ух! , отримаємо (при В-ортогональності векторів 
К-і 
2Ю! та деяких припущеннях щодо матриці демпфування С) систему, яка розпадається 
на незалежні відносно у) рівняння: 


т ' 2 
УРО) я 260 Ук (1) нор УК(І) є РІ), 210, 


) ) 


; (2) 
з уко) 2 У з» 


0 
уко) з У 
-0,5 яТ ротою УС БУТОН о 
де шал", 0«8х1, Ріпузр ВИ), уро зх" Ваг, ур'зх РезАНОТ зді ОР ЛАЕНОНЬ В 
Розв'язки задач (12) мають вигляд: 


ко У? р ва (0) бро (от) | 


Я п в 1 ра 
ук) се"? зітодуі ур собі | з Ї, Р.) зіпа, (1 - тат, 
я 


У 


де 0; - 94/1- 6 . Аналіз цих розв'язків свідчить про те, що істотний внесок в роз- 


в'язок х(7) задачі (8) дають лише близько 10 - 20 складових, які відповідають 
мінімальним власним значенням. 

Існує ряд випадків навантажень, коли можливе точне обчислення розв'язків 
задач (12), наприклад, для вітрового, сейсмічного або гармонійного навантажень. 
В решті випадків розв'язки у (7) обчислюються чисельно, наприклад, методом скін- 
ченних різниць за схемою Ньюмарка. Також задачі Коші (12) можуть розв'язуватися 
методом Рунге - Кутта четвертого порядку, орієнтованим на розв'язування систем 
звичайних диференціальних рівнянь другого порядку. 

Розв'язування узагальненої алгебраїчної проблеми власних значень. В задачі 
на власні коливання, а також, як зазначалось вище, при використанні методу розви- 
нення за формами власних коливань, як правило, необхідно знайти лише невелику 
кількість власних форм (порівняно з порядком задачі), які відповідають мінімальним 
власним значенням. Тому відповідна часткова проблема власних значень (10) може 
розв'язуватися методом ітерацій на підпросторі, паралельні алгоритми якого опи- 
сано в |31. 

Для отримання наближеного розв'язку за допомогою цього методу важливо 
оптимізувати вибір розмірності підпростору, що ітерується, залежно від кількості 
необхідних власних пар. Важливо також з огляду на те, що матриця маси В задачі (10) 
часто є напівозначеною, так вибирати початковий підпростір, щоб матриці проекцій 
на підпростір були додатно означеними. 

Слід також брати до уваги, що кількість ненульових елементів матриці мас В 
задачі (10), як правило, значно менша за кількість елементів матриці жорсткості А, 
які беруть участь в обчисленнях. Тому з точки зору підвищення ефективності 
паралельного алгоритму доцільно зберігати всі ненульові елементи матриці мас, не 
зважаючи на те, що вона є симетричною. Такий підхід значно спрощує паралельний 
алгоритм множення розрідженої симетричної матриці на прямокутну матрицю 1 
зменшує кількість обмінів даними та обсяг цих даних. 

При виборі паралельного алгоритму необхідно враховувати його ефективність. 
Коефіцієнти прискорення та ефективності паралельних алгоритмів, які згадано вище 
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та які можуть використовуватись в задачах розрахунку міцності будівельних кон- 
струкцій, подано в монографії |3| та роботі |З8|. 

На вибір паралельних алгоритмів при розв'язуванні задач з розрідженими мат- 
рицями суттєво впливає також збалансованість завантаження паралельних процесів 
обчисленнями. Тому для кожної задачі має обиратися такий метод оптимізації струк- 
тури матриць 1 такий алгоритм розв'язування задачі та його параметри, які забезпе- 
чують не тільки зменшення загальної кількості арифметичних операцій, але 1 мінімі- 
зацію часу очікування даних окремими процесами при обчисленнях. 


Достовірність наближених розв'язків 


При розрахунку міцності конструкції важливо гарантувати достовірність отри- 
муваних розв'язків. Коли йдеться про достовірність розв'язку, то мається на увазі, 
що досліджується достовірність наближених розв'язків математичних задач (1) - (3). 
Повна похибка наближеного розв'язку враховує вплив похибок вихідних даних, по- 
хибок дискретизації і формування даних відповідної дискретної задачі (8), (9) або (10) 
та похибок комп'ютерного розв'язку цієї дискретної задачі. 

Достовірність розв'язку прикладних задач з наближеними вихідними даними 
гарантується використанням стійких до збурення вихідних даних математичних 
моделей конструкцій, що розраховуються. Використання при дискретизації МСЕ 
теоретично обгрунтованих скінченних елементів, які задовольняють умовам збіж- 
ності і для яких отримано оцінки похибок розв'язку, є передумовою отримання досто- 
вірного наближеного розв'язку задачі. Для отримання достовірного розв'язку також 
важливо узгоджувати точність квадратурних формул, які використовуються при об- 
численнях (7) елементів матриць (жорсткостей, мас, демпфування) і вектора наван- 
тажень, з точністю МСЕ, враховуючи при цьому поведінку підінтегральних функцій. 

Похибка комп'ютерного розв'язку хХ відносно точного розв 'язку х відповід- 
ної дискретної задачі (8), (9) або (10) залежить від точності завдання вихідних даних, 
методу і алгоритму розв'язування цієї дискретної задачі, а також похибок заокруг- 
лення. Достовірність розв'язку також залежить від властивостей комп'ютерної зада- 
чі, які можуть суттєво відрізнятися від властивостей вихідної задачі. Наприклад, при 
розв'язування СЛАР важливу роль відіграє число обумовленості матриці, а при 
розв'язуванні часткової проблеми власних значень необхідно пересвідчитись, що всі 
обчислені власні значення є найменшими. Тому апріорне та апостеріорне досліджен- 
ня математичних властивостей комп'ютерних версій дискретних задач (8), (9) або 
(10) є необхідною умовою отримання достовірного розв'язку. Комп'ютерні методи 
та паралельні алгоритми дослідження властивостей дискретних задач та їх розв'язків 
описано в монографії |31. 

Після отримання розв'язку задачі необхідно провести апостеріорний аналіз 
результатів з точки зору задоволення тим гіпотезам, які були використані при побу- 
дові математичної моделі задачі, а також вимогам точності розв'язку. Якщо роз- 
в'язок задовольняє всім поставленим вимогам, то можна вважати, що цьому роз- 
в'язку можна довіряти і використовувати. 


Технологія дослідження та розв'язування задач розрахунку міцності 
конструкцій на комп'ютерах з паралельною організацією обчислень 

Основна ідея сучасного підходу - інтеграція в рамках єдиного програмного 
комплексу виконання наступних завдань: 


- інтерактивна постановка задачі - створення скінченноелементної моделі 
конструкції; 
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- автоматичне формування даних дискретної задачі; 
автоматичне дослідження та розв'язування дискретної задачі; 

- автоматична обробка результатів розв'язування дискретної задачі; 

- інтерактивний аналіз та використання результатів розрахунку. 

Паралельні обчислення варто використовувати для завдань, які виконуються 
автоматично, в першу чергу, для дослідження та розв'язування дискретних задач 
(тому що, як правило, це завдання вимагає найбільше часу). 

Як наголошувалося вище для дослідження та розв'язування задач (1) - (3) на 
комп'ютерах з паралельною організацією обчислень створено програмний комплекс 
(ПК) розв'язування задач розрахунку міцності будівельних конструкцій |21, 31. Цей 
ПК (надалі ПК ЛІРА-сСІизіег) розроблено на базі ПК ЛІРА та інтелектуальних пара- 
лельних програм бібліотеки Ірагіїб |2), |3| для інтелектуальних паралельних ком- 
п'ютерів Інпарком. Після деякої адаптації (яка в більшості випадків полягає у пере- 
компіляції паралельних програм) ПК ЛІРА -сіцзісг може бути встановлено на інших 
багатопроцесорних комп'ютерних системах. 

При створенні цього програмного комплексу необхідно було взяти до уваги, що 
графічний інтерфейс пре- та постпроцесора ПК ЛІРА, як і в більшості аналогічних 
програмних комплексів, працює в середовищі ОС У/тпдоууз, а програмні модулі, які 
виконують паралельні обчислення, - в ОС Ілпих (ця ОС використовується на біль- 
шості багатопроцесорних комп'ютерів). Тому при створенні комплексу використано 
те, що в склад багатопроцесорної робочої станції Інпарком входить По8і-система, на 
одному з комп'ютерів (графічній станції) якої може бути завантажено ОС У/іпдомуз, 
в той час як на обчислювальні вузли - ОС Ілпих |2|. Як ро8і-комп'ютер в ПК ЛІРА- 
сПи5їег може використовуватись комп'ютер користувача, зокрема персональний. Та- 
кий підхід дозволяє поетапно переходити до використання паралельних обчислень 
при розрахунку міцності будівельних конструкцій, починаючи з найбільш ресурсо- 
ємного етапу - дослідження та розв'язування дискретних задач. 

Розрахунок міцності конструкції за допомогою ПК ЛІРА-сіизієг проводиться 
наступним чином: 

- на графічній станції Інпарком інтерактивно на мові предметної області, 
використовуючи препроцесор ПК ЛІРА, створюється або вводиться раніше створена 
скінченноелементна модель конструкції; 

- на графічній станції автоматично препроцесором програмного комплексу 
формуються вихідні дані для формування (можливо, часткового) та розв'язування на 
обчислювальних вузлах відповідної дискретної задачі (3) або (10); ці дані запи- 
суються в бінарні файли, використовуючи прийняті в ПК ЛІРА формати даних; 

- автоматично за допомогою спеціального інтерфейсу на обчислювальні вузли 
комп'ютера Інпарком завантажується паралельна програма, яка: 1) визначає алго- 
ритм і параметри розв'язування задачі, 2) зчитує вихідні дані з файлів та формує 
розподілену між паралельними процесами відповідно до вибраного алгоритму роз- 
в'язування дискретну задачу, 3) досліджує і розв'язує дискретну задачу, викори- 
стовуючи відповідні модулі бібліотеки Шраті, 4) зберігає результати дослідження 1 
розв'язування задачі у текстовому та бінарних файлах у форматах ПК ЛІРА для по- 
дальшої обробки і використання постпроцесором програмного комплексу; 

- на графічній станції зчитуються з файлів результати розв'язування дискрет- 
ної задачі та автоматично засобами ПК ЛІРА проводиться обробка отриманих ре- 
зультатів; 

- на графічній станції інтерактивно на мові предметної області постпроцесором 
ПК ЛІРА проводиться аналіз та використання результатів розрахунку або результати 
зберігаються для подальшої обробки та використання на комп'ютері користувача. 
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Для розв'язування динамічних задач використовується метод розвинення за 
формами власних коливань. При цьому дослідження та розв'язування часткової уза- 
гальненої алгебраїчної проблеми власних значень виконується паралельно на об- 
числювальних вузлах, а розв'язування рівнянь (12) реалізується ПК ЛІРА, для чого в 
бінарних файлах потрібно зберігати обчислені власні значення та вектори, а також 
трикутне розвинення матриці жорсткостей. 

У випадку розв'язування нелінійної задачі ітераційний процес побудовано так, 
що на кожній ітерації на обчислювальних вузлах для чергового наближення форму- 
ється матриця жорсткостей задачі (8) та виконується Її трикутне розвинення, а всі 
інші дії здійснюються ПК ЛІРА. Тобто для отримання розв'язку необхідно багато- 
разово звертатися до обчислювальних вузлів багатопроцесорного комп'ютера. В по- 
дальших версіях ПК ЛІРА-сіцзіег планується повністю організувати ітераційний 
процес на обчислювальних вузлах, але це вимагає внесення значних змін у ПК ЛІРА. 

Інтерфейс користувача ПК ЛІРА-сішцяіег. В ПК ЛІРА-сіцзісг використову- 
ється інтерфейс ПК ЛІРА, який доповнено двома вікнами. Тому робота користувача 
ПК ЛІРА -сіизіег практично не відрізняється від роботи з ПК ЛІРА. 

Перше з цих вікон - стартове вікно ПК ЛІРА-сіизіег (рис. 1), з відкриття якого 
розпочинається робота з програмним комплексом. У ньому міститься деяка довідко- 
ва інформація про встановлений режим роботи програмного комплексу - кластерний 
(з паралельним процесором) або стандартний (з послідовним процесором) - та про 
стан багатопроцесорного комп'ютера (кластера) - дані про загальну та вільну 
кількість ядер. Також тут містяться кнопки, що дозволяють розпочати або закінчити 
сеанс роботи з програмним комплексом, змінити режим його роботи, відновити 
вікно вибору (зміни) кількості ядер (рис. 1). 


в) ЛИРА-Физіег 
Информация Состояние кластера 


Добро пожаловать в программньій комплекс са дано С 


ЛИРА-сіийзіег для работь с ПК ЛИРА вер.9.4, Свободнье ядра ГР 
используя интеллектуальную рабочую 


станцию Инпарком (кластер). 3 
бновить | 


| Установлен стандартньй режим ПК ЛИРА 


Перейти в кластерньїй режим ПК ЛИРА Начать сеанс ЛИРА-сійзівг 
Завершить сеанс ЛИРА»-сіцзівг 


Рисунок І - Стартове вікно ПК ЛІРА-сІизіег 


Друге вікно - вікно зміни кількості ядер - пов'язано безпосередньо з самим 
процесом паралельних обчислень. Так, коли закінчується формування даних для 
роботи паралельної частини ПК ЛІРА-сІи5зі»г, відкривається це вікно. У ньому вка- 
зується загальний ресурс та кількість вільних ядер, яка визначається відповіддю на 
спеціальний запит кластеру безпосередньо відкриттям цього вікна. Також користу- 
вачу пропонується (поле «Рекомендовано») автоматично визначена, найкраща для 
даної задачі (з точки зору швидкодії) кількість ядер. У випадку, якщо користувача не 
задовольняє рекомендована кількість, він власноруч може задати інше число на свій 
розсуд. Введене число перевіряється на коректність, у тому числі й з урахуванням 
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реальних потужностей кластера. Для зручності у вікні розміщено кнопку оновлення 
кількості процесів. Якщо на даний момент кількість вільних ядер замала для даної 
задачі, користувач може очікувати на більше число, час від часу натискаючи кнопку 
оновлення. Коли користувач вважає, що задачу можна запускати на кластері, він 
натискає на кнопку «ОК», вікно зникає і розпочинається підготовка до запуску 
паралельної програми. Якщо ж з моменту появи вікна користувач не робив спроб 
змінити кількість ядер, через деякий час (інформація про який також дається у вікні) 
вікно закривається і розрахунок розпочинається автоматично. 

Користувач також має можливість відмовитись від подальшої появи вікна ви- 
бору кількості ядер. В такому випадку це число визначається автоматично. Для від- 
новлення режиму появи цього вікна необхідно використати відповідну кнопку стар- 
тового вікна ПК ЛІРА-сіІи5іег (яка активізується у момент відмови). Вікно вибору 
кількості ядер містить також поле для підказок, які подаються при наведенні курсора 
миші на відповідні елементи вікна. 


ні Изменение количест... | соні! то 


Всего ядер 258 
Рекомендовано 0 


Свободнье ядра 8 


Использовать 0 


3 Обновить 


Окно будет закрьто через 3 секунд 


Г' Больше не показьквать окно в зтом сеансе 


Количество ядер, которье будут 
использовани для решения задачи на 
кластере 


Рисунок 2 - Вікно вибору (зміни) кількості ядер 


Перед запуском відповідної паралельної програми проводиться редагування 
команди запуску - встановлюється кількість процесів (яка дорівнює вибраній 
кількості ядер), якщо це число відмінне від наявного в команді запуску. 

Користувач може працювати з ПК ЛІРА-сіІц5іег на інтелектуальному паралель- 
ному комп'ютері Інпарком - безпосередньо або використовуючи підключення через 
Інтернет до віддаленого робочого столу. Також передбачена можливість роботи з 
програмним комплексом, використовуючи як Бо8і-комп'ютер персональний комп'ю- 
тер користувача та Інтернет для зв'язку з багатопроцесорним комп'ютером. В остан- 
ньому випадку перед обмінами даними між вказаними комп'ютерами проводиться їх 
архівація з метою зменшення обсягу даних, які пересилаються засобами каналів 
зв'язку. 

Програмним комплексом ЛІРА-сіц5їег на робочих станціях Інпарком прово- 
дилося розв'язування тестових задач. При цьому час розв'язування системи лінійних 
алгебраїчних рівнянь (8) зменшувався в 10 - 20 раз, а час розв'язування часткової 
алгебраїчної проблеми власних значень (10) - в 5-10 раз порівняно з розв'я- 
зуванням на персональному комп'ютері. Деякі приклади наведено в роботах (5, 19). 
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Висновки 


Проведені авторами дослідження і апробація програмного комплексу для роз- 
в'язування задач розрахунку міцності будівельних конструкцій на робочих станціях 
Інпарком показали перспективність створення програмних засобів для розв'язування 
прикладних задач на паралельних комп'ютерах кластерного типу шляхом вбудо- 
вування паралельних програм для розв'язування окремих підзадач в наявні прог- 
рамні комплекси в різних галузях науки 1 техніки. Як такі паралельні програми мо- 
жуть використовуватися програми бібліотеки Прагії або її аналогів на інших ком- 
п'ютерах з паралельною обробкою інформації. 

Також є перспективним розширення кола задач, які розв'язуються за допо- 
могою програмного комплексу ЛІРА -сІизіег на робочих станціях Інпарком, зокрема, 
більш глибокого розпаралелювання розв'язування нелінійних задач, а також обробки 
результатів розв'язування дискретних задач. Заслуговує на увагу також перенесення 
на обчислювальні вузли формування матриць (жорсткостей, мас) і векторів наванта- 
ження окремих скінченних елементів. 
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А.Н. Химич, А.В. Попов, В.В. Полянко 

Технологиий вьтсокопроизводительньх вьтчислений для исследования и решения задач 

расчета прочности конструкций 

Рассматриваются технологий вьісокопроизводительньтх вьтічислений для исследования и решения задач 
(статических линейньх и нелинейньх и динамических линейньх) расчета прочности конструкций. 
Описанью параллельньве алгоритмь: формирования, исследования и решения систем разрешающих 
уравнений, получаємьіх при дискретизации методом конечньтх злементов названньтх задач, а также 
реализация предложенньтх технологий в программном комплексе решения задач расчета прочности 
строительньх конструкций. 


А. Кнітісі, А. У. Ророу, У.У. Роіуапко 

Тесіпоїобіезя ої Нієр-регіогтапсе Соприбабіоп5 їог бе Пуе5біданіпо апа 5о0Їуіпо ої Ргобіет5 

оп Ше 5ігепої Апаїузбіз ої Соп5ігисбіоп5 

Тесфпоїоєіе8 ої Біє|-регіогтапсе согариїайнопзя їог іпуезівайтя апа 50Імїпо, ої ргобіетя (5кайс Ппеаг апа 
поп-Ппеаг, дупатіса! Ппеаг) оп Де 5ітепеЮ апаіубі8 ої соп5ігисійопя аге деаїї мії. РагаПе! аїєогіїпта5 Гог 
Фе Богпаїпє, іпуе5іїсайття апа 50Їміпо, ої 5убієта ої ге5оЇміпе, едцайопя обіаїпед дигіпе їбе Ппіе-еЇетепі 
Фізсгейігайоп ої Фе П5кед ргобієтія а5 ууе/! а5 ппріетепіайоп ої іесрпоїогіе5 ргорозед іп Ше ргобгат 
сотріех Гог Фе 50ЇІміпо, ої ргобіета8 оп Бе з5ігепеїЮ апаіузіз ої соп5ігисіоп5 аге Фезстібед. 
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